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The electronic structure of an atom with Z < Z c = 137 can be described by the Dirac equation 
with the Coulomb field of a point charge Ze. It was believed that the Dirac equation with Z > Z c 
meets difficulties because the formula for the lower energy level of the Dirac Hamiltonian formally 
gives imaginary eigenvalues. But a strict mathematical consideration shows that difficulties with the 
electronic spectrum for Z > Z c do not arise if the Dirac Hamiltonian is correctly defined as a self- 
adjoint operator. In this article, we briefly summarize the main physical results of that consideration 
in a form suitable for physicists with some additional new details and numerical calculations of the 
electronic spectra. 
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INTRODUCTION 



The question of electronic structure of an atom with large charge number Z of the nucleus, especially with Z that is 
more than the critical value Z c — oT 1 ~ 137, 04, where a is the finite structure constant, is of fundamental importance. 
The formulation of QED cannot be considered really completed until an exhaustive answer to this question is given. 
Although nuclei with overcritical charges can hardly be synthesized (at present, the maximum is Z = 118), the 
existing heavy nuclei can imitate the supercritical Coulomb fields at collision. Nuclear forces can hold the colliding 
nuclei together for 10~ 19 s or more. This time is enough to effectively reproduce the experimental situation where the 
electron experiences the supercritical Coulomb field [2|. 

The electronic structure of an atom with Z < Zr can be described by the Dirac equation, which gives relativistic 
■ electronic spectra in agreement with experiment [3j. For such Z a complete set of solutions of the Dirac equation 
J> 1 exists; the corresponding Coulomb field does not violate the vacuum stability, therefore, the Furry picture in QED 
can be constructed, and the relativistic quantum mechanics of an electron in such a Coulomb field based on the Dirac 
equation is self-consistent. It was believed that the Dirac equation with Z > Z c meets difficulties 041] • One of the 
\ standard arguments is that the formula for the lower energy level, 

E = mc 2 yjl- (Zaf , (1) 

formally gives imaginary result for Z > Z c . This difficulty of the imaginary spectrum was attributed to an inadmissible 
singularity of the supercritical Coulomb field at the origin for a relativistic electron, see Q • It was believed that this 
difficulty can be eliminated if a nucleus of some finite radius R is considered. It was shown that with cutting off the 
Coulomb potential with Z < 173 at a radius i?~ 1,2 x 10~ 12 cm, the Dirac Hamiltonian has physically meaningful 
spectrum and eigenstates But even in the presence of the cutoff, another difficulty arises at Z ~ 173. Namely, the 

lower bound state energy descends to the upper boundary E = —mc 2 of the lower continuum, and it is generally agreed 
that in such a situation, the problem can no longer be considered a one-particle one because of the electron-positron 
pair production, which, in particular, results in a screening of the Coulomb potential of the nucleus. Probabilities 
of the particle production in the heavy- ion collisions were calculated within the framework of this conception Q. 
Unfortunately, experimental conditions for verifying the corresponding predictions are unavailable at present. 

Not disputing the fact that taking account of a finite size of the nucleus corresponds to a more realistic setting up 
the problem, we do not agree with the assertion that the Dirac Hamiltonian with the Coulomb field of overcritical 
point-like nucleus charge is inconsistent. The above-mentioned difficulties with the spectrum for Z > Z c do not 
arise if the Dirac Hamiltonian is correctly defined as a self-adjoint (s.a.) operator. A first heuristic attempt in this 
direction is due to 0. A rigorous mathematical treatment of all the aspects of this problem including a spectral 
analysis of the Hamiltonian based on the theory of s.a. extensions of symmetric operators and the Krein method 
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of guiding functionals was presented in [1, was demonstrated that from a mathematical standpoint, a 
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definition of the Dirac Hamiltonian as a s.a. operator presents no problem for arbitrary Z. A specific feature of the 
overcritical charges is a nonuniqueness of the s.a. Dirac Hamiltonian, but this nonuniqueness is characteristic even for 
Z > Z s = (y/3/2) a -1 ~ 118, 68. For each Z > Z s , there exist a family of s.a. Dirac Hamiltonians parametrized by a 
finite number of extra parameters (and specified by additional boundary conditions at the origin). The existence of 
these parameters is a manifestation of a nontrivial physics inside the nucleus. A real spectrum and a complete set of 
eigenstates can be evaluated for each Hamiltonian, so that a relativistic quantum mechanics for an electron in such a 
Coulomb field can be constructed. 

In the present article, we briefly summarize all the previously obtained formal results in a form more suitable 
for physicists with some additional new important details and numerical calculations of the electronic spectra. The 
spectrum of each Hamiltonian consists of a universal continuous part that is a union of two intervals (— 00, m] and 
[m, 00) and a specific discrete spectrum located in the interval [— m, m). We concentrate on the discrete spectrum. 



2. DIRAC HAMILTONIAN WITH COULOMB FIELD 



We consider an electron of charge — e < and mass m moving in the Coulomb field of charge Ze > 0. We describe 
this field by a scalar electromagnetic potential of the form Aq — Zer^ 1 , we set h = c= 1 in what follows. A behavior 
of the electron in the Coulomb field is governed by the Dirac Hamiltonian H (Z) that is a s.a. operator in the Hilbert 
space fj of square-integrable bispinors \& (r) . On its domain which must be properly specified, H (Z) acts by the 
differential operation 



H (Z) = 7 (7P + m) — qr 1 , p = — iV, r = |r| , q = Za 



(2) 



(in what follows, we use 7-matrices in the standard representation). In the problem under consideration, there are 

-2 - 
three commuting s.a. operators J , J z , and K, where J is the total angular momentum and K is the so-called spin 

operator 



J = L + S/2, L = [r x p] , K = 7 1 + (sl) 



All they commute with H (Z). Any bispinor <]/ (r) can be represented as \l/(r) = J^j m c ^jMX ( r )' where ^fj,M,c are 
bispinors of the form 



(3) 



&j,M,C are normalized spherical spinors, / (r) and g (r) are radial functions, and j = 1/2, 3/2, M — —j, — j + j, 
C = ±. Bispinors ^j,M,C, are eigenvectors of J , J z , and K, 
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J$ = j(j + 1)*, J z ^ = M*, K'fy = -CO' + 1/2)* . 

Let L 2 = L 2 (K + ) L 2 (IR + ) (where L 2 (R + ) is the space of functions of r square-integrable on the semiaxis 
4. = [0,oo)) be the Hilbert space of doublets F(r), 



F(r) 



( m \ 



(f(r)/g(r)) 



with the scalar product 

(F U F 2 ) 

Then rep. §3§ and the relation 



drF+ (r) F 2 (r) 



dr 



fi (r)/2 (r) +gi (r)g 2 (r) 



I* 



dr [|/(r)| 2 + | 5 (r)| 5 



show that any subspace of bispinors *j,M,i with fixed j, M, ^ is unitary equivalent to L 2 (K+), an explicit form of this 
equivalence is 



*j,Af,c(r) = 7- _1 n j ,M,c(e,¥')F(r), F(r) = r J sin 6McbpUf MX (9, ip)Vj, M , C (r). 



(4) 
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Here, Tlj.MX and ^ are the respective (4 x 2)- and (2 x 4)-matrices, 



J sin BMd/p [nt M iC (0, <^)n J - MiC (0, 



= ^afc, a, 6 = 1,2, 

26 



where = (0/0) is a two-column and T = (0, 0) is a two-row. The stationary Schrodinger equation H (Z) \fr (r) = 
Ety (r) in the Hilbert space F) is reduced to radial equations 

h(ZJ,t)F(r) = EF(r), feL 2 (i + ), 

in the Hilbert space L 2 (R + ), where h(Z,j,Q are s.a. partial radial Hamiltonians acting on the doublets F (r) by 
radial differential operations 

h(Z,j,() = -i(T2d r + C(j + l/2)r- 1 a 1 - qr' 1 +m<r 3 . (5) 

The problem of constructing a rotationally invariant s.a. Dirac Hamiltonian H (Z) is reduced to the problem of 
constructing s.a. partial radial Hamiltonians h(Z,j,Q- We construct all possible s.a. partial radial Hamiltonians 
h(Z,j,£) using the theory of s.a. extensions of symmetric operators, which is reduced to specifying their domains 
of definition Du% j q C L 2 (R + ). By construction, each operator h(Z,jX) is a s.a. extension of the so-called initial 
symmetric operator h m (Z,j, £) with the domain , z ■ q — 2?(R + ) 2?(R + ), where 2?(R+) is a space of smooth 
compactly supported functions on the semiaxis R+, and (simultaneously) is, generally, a s.a. restriction of the 
adjoint operator kt (Z,j, Q defined on the so-called natural domains Di . (R+) of doublets F G L 2 (R + ) that are 

absolutely continuos in R + and such that h (Z,j, C) F € L 2 (R + ). Thus, D~ h ^ z . ^ C D^ z ^. ^ C . ^ (R+), and 

if j () does not coincide with D* h ^ z . ^ (R+), it is specified by some additional asymptotic boundary conditions 
at the origin, which are defined not uniquely. 

A result of constructing s.a. radial Hamiltonians h(Z,j,() essentially depends on the values of the parameters Z 
and j. 

There are two regions in the first quadrant of the j, Z plane, we call them the nonsingular and singular ones, where 
the problem of s.a. extensions has principally different solutions. These regions are separated by the singular curve 
Z = Z s (j), where 



z s (j) = vO'0' + i) 



a . 



so that the nonsingular and singular regions are defined by the respective inequalities Z < Z s (j) and Z > Z s (J). The 
values 

Z s (1/2) = 118,68; Z s (3/2) = 265, 37; Z s (5/2) = 405, 36; 

can be called the singular values of Z for a given j, see FIG. [TJ 

In what follows, we describe s.a. radial Hamiltonians h(Z,j,^) and their spectra in the nonsingular and singular 
regions separately. The consideration in the non-singular and singular regions is true for the whole set of all quantum 
numbers. 

All the mathematical details concerning s.a. extensions and a spectral analysis of Hamiltonians (including complete 
sets of eigenfunctions) based on the Krein method of guiding functionals can be found in [l|, [ll|, [l2[ • 



NONSINGULAR REGION 



In the nonsingular region, Z < Z s (j), the partial radial Hamiltonian h{Z,j,C l ) is defined uniquely, h(Z,j,£) 
(Z,j, £), and its domain is j q — ^h(z j f) 0^+)' ^he functions belonging to D q (®+) have the following 



asymptotic behavior 

F(r) = 0(r 1/2 ), r ->■ 0; F{r) -> 0, 
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FIG. 1: Regions in the j, Z plane. Singular and critical curves. 

A discrete spectrum {E n (o (Z, j, £)} of each Hamiltonian h (Z, j, £) is given by 

E n(0 (Z, j, = ^(" <C) +7) = v / C? + i/ 2 ) 2 - (72 . (6) 

vV + (n(0 + 7) 2 v 

The quantum numbers take the values 

G N = {1, 2, ...} ; G Z+ = {0, 1, 2, ...} . (7) 

The expression ([6]) coincides with the well-known Sommerfeld formula for the spectrum of the Dirac electron in the 
Coulomb field. This result justifies the standard formal treatment of the Dirac Hamiltonian with Z in the nonsingular 
region in the physical literature where the Dirac Hamiltonian is identified with the differential operation Hd(Z) ([2]) 
and the natural domain is implicitly assumed. 



4. SINGULAR REGION 

In the singular region, Z > Z s (j), the s.a. radial Hamiltonian h(Z,j,() with given Z,j,( is defined not uniquely 
(the deficiency indices of the operator hi a (Z,j, () are (1, 1)). There exists a one-parameter family{/i„ (Z,j, £) , v € 
[ — 7r/2, 7r/2], — tt/2 ~ tt/2} of Hamiltonians denned on different domains D~ h (z j *- ^h{q «•) 0^+) s P ec i ne d by 
different asymptotic s.a. boundary conditions at the origin, which are parametrized by the parameter v. 

The position of the discrete energy levels E^ () {Z, j, £) essentially depends on is, in particular, there exists a value 
v = is- rn , for which the lower energy level coincides with the boundary E = —m of the lower continuous spectrum. 

Technically, it is convenient to divide the singular region into three subregions, we call them subcritical, critical, 
and overcritical regions. The subregions are distinguished by a character of asymptotic boundary conditions at the 
origin. 

The boundary conditions are similar in each subregion, which provides similar solutions of the corresponding spectral 
problems. In what follows, we describe these subregions, the domains D~ h ^ in these subregions, and some details 
of discrete spectra. 



4.1. Subcritical regions 



The subcritical region is defined by the inequalities Z s (j) < Z < Z c (j) , where 

Z c (j) = (j + 1/2) a' 1 , 



(8) 
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see FIG. [TJ The values 

Z c (1/2) = Z c = 137, 04; Z c (3/2) = 274, 08; Z c (5/2) = 411, 12; ... 

can be called the critical Z- value for a given j. In the subcritical region, the s.a. radial Hamiltonians h v (Z,j, £) are 
specified by s.a. boundary conditions, 

F(r) = c[(mr) 7 d+ cos v + (mr)" 7 (i_ sin v] + 0(r 1/2 ), r -> 0, (9) 



where 7 = y (j + 1/2) 2 - <? 2 , < 7 < 1/2, c is an arbitrary complex number, and d± are some constant doublets. 

The spectrum of each h u (Z, j, £) is simple (nondegenerate) and consists of a continuous part that is the set 
(— 00, — to] U [to, 00) and a discrete part located in the interval [—to, to). The discrete spectrum is a growing infi- 
nite sequence {E^ (Z,j, £)} of the energy levels that are the roots of the equation 

/(£)cos^ + r(l-2 7 )sin^ 

/(£)sin*/-r(l-27)cos!/ ' [ ' 

r(l + 2 7 )r(- 7 - qEr^Mm - E) - (x + 7 )r] 
M j r(7- gJ B/r)[g(m- J E)-(x-7)r](2r/TO) 27 ' 

» = CO' + 1/2), r = Vm 2 -E 2 , 

the integers take the previous values ([7]). 

We outline the most important features of the discrete spectrum. For each Z,j, (, there exists v = v~ m = v_ m (Z,j) 
(independent of () such that the lowest energy level is equal to —to. This v- m is determined from eq. (fT0|) by setting 
E=-m and noting that f(-m) = T(l + 2-y){2q)- 2 ~< to yield 

tan ,_ m = _£(i±^) (2(Z )- 2 7 
r(i-2 7 ) v q; 

For fixed E ? () (Z,j,Q as functions of v are monotonically decreasing functions with the properties 



iy-s--7r/2+0 

Let 



lim £7^. (Z, j,C) = Um J5 l 2, , . (Z, j,0 = £„c 



to, £ (o ), there are no energy level for v G (i/_ m ,7r/2], while for each ^ e [— 7r/2, u- m ], there is one level 



_(0 _ n (0 _ J !. C = + 

"0 — "min — [ £ = — ' 

A subtlety is that the function E^ C) (Z,j,Q is defined only for v 6 [— 7r/2, V— m ], which implies that in the energy 
interval 

E^Q (Z, j, C) monotonically growing from —to to £ <o when v changes from v- m to — tt/2. 

The functions E^ {Z 1 j,C t ) with n^) > are defined for all v € [— 7r/2, tt/2]. In each energy interval 

[£ n (0 7 £ n (o+i] , there is one level (Z, j, £) monotonically growing from £ n (o to £„(o+i when changes from 7r/2 

to — 7r/2. We note that the states with the energies E^' 2 ' (Z,j,Q and E^ 2 ^ (Z,j,() — E^/ 2 ^ (Z,j,() = £ n (o 
represent the same eigenstate. It follows from the fact that according to eq. ([9]), the values v — tt/2 and v = —tt/2 
are equivalent and therefore the Hamiltonians hw2 (Z,j,£) and /i_ w /2 {Z,j,Q are the same. The eigenvalues £„<o 
can be found explicitly. For v — ±ir/2, eq. (fT0|) is reduced to the equation 1// (E) — 0, and we find 

(n^' — 7)771 



In particular, we see that the discrete spectrum E ^ {Z,j, Q is given by eq. (jSJ) with Z corresponding to the region 
under consideration. 

For illustration, we give graphs of five low energy levels [Z = 121, j = 1/2) as functions of v, for C, = + (FIG. Oi), 
for ( = — (FIG. [5b), and also a graph of the parameter V- m \j=i/2 as a function of Z, see FIG. [3b. 
In addition, to give an idea of energy scale, we give Table 1 of numerical values of some energy levels. 



FIG. 2: ^-dependence of energy levels E^ v ) a (121, 1/2, £ = ±) and Z-dependence of v- 



\j=i/2 ■ 



Table 1. Energy levels for Z = 121, j = 1/2 



1/ 


Eo 


Ei 


E 2 


E3 


E 4 


"?, 




-0.469411 


0.515067 


0.515067 


0.866198 


0.866198 


0.944175 


0.944175 


0.970121 


0.970121 






-0.487899 


0.465810 


0.514264 


0.865181 


0.866078 


0.943927 


0.944142 


0.970026 


0.970108 











0.46941 


0.85715 


0.85715 


0.941615 


0.941615 


0.969107 


0.941615 
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-0.450991 


0.563241 


0.515898 


0.867497 


0.866322 


0.944482 


0.944209 


2 








-0.469411 


0.515067 


0.515067 


0.866198 


0.866198 


0.944175 


0.944175 


c 


+ 




+ 




+ 




+ 




+ 





4.2. Critical region 



The critical region is the critical curve Z = Z c (j) , see FIG. [TJ For integer Z, this region does not exist if the fine 
structure constant a is an irrational number, see ((HJ. In particular, this region certainly is absent for j = 1/2. In 
the critical region, the s.a. radial Hamiltonian h v (Z,j, £) is specified by s.a. boundary conditions at the origin of the 
form 



F(r) = c [d (r) cosv + d + sinu] + ©(r 1 / 2 lnr), r ->• 0, 



(11) 



where do(r) are some doublet with the asymptotic behavior do(r) = O (In mr) as r — > 0. 

The spectrum of each h v (Z, j, £) is simple (nondegenerate) and consists of a continuous part that is the set 
(—00,— to] U [m, 00) and a discrete part located in the interval [—to, to). The discrete spectrum is a growing infi- 



nite sequence {E X (Z,j, £)} of energy levels that are the roots of the equation 



g(E) cos v — sin v 
g(E) sin v + cos v 



0, g(E) = ln(2r/m) 



+ (j + 1/2) E/r) + 



C - (t + Cm) /E 
2 {j + 1/2) 



-2^(1), 



(12) 



where ip(x) = V (x)/T(x). 

We outline the most important features of the discrete spectrum. For each Z,j, £, there exists v = v- m = v- m (Z,j) 
(independent of Q) such that the lowest energy level is equal to —to. This v^ m is determined from eq. (|12[) by setting 
E = —to and noting that g(—m) = ln(2g c ) — 2?/>(l) + (/q c to yield 

tani/_ m = ln(2g c ) - 2^(1) + (/q c - 
For fixed E (Z,j, Q as functions of v are monotonically decreasing functions with the properties 



lim E^ () (Z,j,0 



lim E%{. )+1 (Z,j,£)=£ nM . 



j/->7r/2-0 
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A subtlety is that the function E^L (Z,j, C) is defined only for v G [— 7r/2, f_ m ], which implies that in the energy 



interval 



level E^ C) (Z,j,Q monotonically growing from — m to £ «) when v changes from v- m to —ir/2. The functions 

n n o 

E^ (Z, j, C) with > rig are defined for all v G [— tt/2, 7r/2]. In each energy interval [£„<o , £ n (0+i] > there is 
one level i?^, (Z,j,() monotonically growing from £ n (o to £„(o+i when v changes from 7r/2 to — tt/2. We note 



-m,£ a) 



, there are no energy level for v G (^_ m ,7r/2], while for each f G [— 7r/2, there is one 



that the states with the energies -E„~o (^>i>0 an d %J<)+i (^,3,0 = -^„7o J, C) = £n«) represent the same 

eigenstate. It follows from the fact that according to eq. (ITT|) . the values f = 7r/2 and v = —tt/2 are equivalent 
and therefore the Hamiltonians h^/ 2 (Z, j, £) and {Z, j, Q are the same. The eigenvalues can be found 

explicitly. For v — ±tt/2, eq. (fl2|) is reduced to the equation 1/g (E) — 0, and we find 



,.(0 



v /(j- + l/2) 2 + (n«)) 1 



(13) 



For illustration, we give graphs of five low energy levels (j = 1/2) as functions of v, for £ = + (FIG. Eh,), for £ = 
(FIG. [5}d). In addition, to give an idea of energy scale, we give Table 2 of numerical values of some energy levels. 
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FIG. 3: f-dependence of energy levels E^ () (Z c , 1/2, ( — ±). 



Table 2. Energy levels for Z c , j = 1/2 





E 


Ei 


E 2 


E 3 


E 4 


"2 







0.707107 


0.707107 


0.894427 


0.894427 


0.948683 


0.948683 


0.970143 


0.970143 


"T 




0.476777 


0.548043 


0.581793 


0.864486 


0.939071 


0.93892 


0.965944 


0.96591 


0.978405 







-0.636533 


0.466006 


0.539248 


0.854923 


0.856797 


0.936069 


0.936298 


0.964727 


0.964778 


"t 






0.290246 


0.476777 


0.843971 


0.803545 


0.932388 


0.920568 


0.96309 


0.963161 


1 











0.707107 


0.707107 


0.894427 


0.894427 


0.948683 


0.948683 


C 


+ 




+ 




+ 




+ 




+ 





4.3. Overcritical region 

The overcritical region is defined by the inequality Z > Z c (j), Z c (j) = (j + 1/2) or 1 . 

In this region, the s.a. radial Hamiltonians h v (Z,j,£) are specified by s.a. boundary conditions 

F(r) = c [ic'>r)'> + - ie~ iv (mr)'™ p.] + 0(r 1/2 ), r -> 0, 



(14) 



where a = y q 2 — (j ' + 1/2) , c is an arbitrary complex number, and p± are some constant doublets. 

The spectrum of each h v (Z, j, £) is simple (nondegenerate) and consists of a continuous part that is the set 
(— oo, — to] U [to,oo) and a discrete part located in the interval [— m, m). The discrete spectrum is is a growing 
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infinite sequence {-E« (Z,j, £)}, n G Z+, of energy levels that are the roots of the equation 



cos[e(£) -v]=Q, 



(15) 



(lnS a ) 



(7 In • 



2r 



where fli = -2ia, B 2 (E) = ia - Egr -1 , and B 3 (£) = r(j + 1/2 - <cr) - C?(m - E). 

We outline the most important features of the discrete spectrum. For each Z,j, £, there exists v = v- m = i>- m (Z,j) 
(independent of Q such that the lowest energy level is equal to — m. This v^ m is determined from eq. (|15[) by setting 
E = —m which yields 



2*e(-m) 



T(-2i(r) 
T(2ia) 



(2?) 



2>cr 



The energy levels are determined by 

e (-m) - e (z,j,c) 



For fixed n, (Z, j, £) as functions of v are monotonically decreasing functions with the properties 



lim EM{Z,j,0 

>-tt/2+0 



}^_ o EiU^x) 



A subtlety is that the function Eq (Z,j, £) is defined only for v 6 [— tt/2, U- m ], which implies that in the energy 



E 7r ^ 2 ' ) (Z,j,Q^, there are no energy level for v 6 (z/_ m , 7r/2], while for each v € [— 7r/2, U- m ], there is 

(i?, j, C) monotonically j 
The functions E^ ] (Z, j, C) with 



interval 

one level Eq^ (Z,j, Q monotonically growing from —m to Eq lr ^ 2 ' 1 (Z,j, £) when v changes from v_ m to —tt/2 



n > 



E^' 2) (Z, j, = E^' 2) (Z, j, , E^ /2) (Z, j, o 



1 are defined for all v € 

there is one level E y n 



E^ /2) (Z,j,0 = E [ n _\ /2) {Z,j,Q to E { n " /2) (Z,j,C) when v changes from tt/2 to 



[— tt/2, tt/2]. In each energy interval 
1 (Z, j, C) monotonically growing from 

7r/2. We note that the states 
with the energies E^^^ 2 ^ [Z,j,Q) and E^/*' (Z,j,Q represent the same eigenstate. It follows from the fact that 
according to eq. (fH|) . the values v — tt/2 and v = —tt/2 are equivalent and therefore the Hamiltonians h w /2 (Z,j, £) 

and h-tii (Z,j, £) are the same. 

We note that in contrast to previous regions, the integers n take the same values for both £ = + and £ = — because 
there is no the Sommerfeld degeneracy. 

For illustration, we give graphs of five low energy levels [Z = 138, j = 1/2) as functions of v, for £ = + (FIG. |4^,), 
for C = — (FIG. 2b), and also a graph of the parameter V-m\j=\/2 as a function of Z, see FIG. 0fc- In addition, to 
give an idea of energy scale, we give Table 3 of numerical values of some energy levels 

Table 3. Energy levels for Z = 138, j = 1/2 





Eq 


Ei 


E2 


E3 


E 4 


~ ~ 


0.705525 


0.024086 


0.893927 


0.717081 


0.948327 


0.89653 


0.969889 


0.949229 


0.980395 


0.970294 


"T 


0.663242 


-0.075165 


0.884039 


0.682943 


0.944912 


0.887947 


0.968366 


0.946203 


0.979595 


0.96893 





0.281658 


-0.42001 


0.839196 


0.584008 


0.931198 


0.865425 


0.931198 


0.938723 


0.976677 


0.965675 


T 






0.742474 


0.165125 


0.903704 


0.764448 


0.951885 


0.909191 


0.971521 


0.95389 


2 






0.705525 


0.024086 


0.893927 


0.717081 


0.948327 


0.89653 


0.969889 


0.949229 


c 


+ 




+ 




+ 




+ 




+ 





Table 4. Energy levels for Z = 180, j = 1/2 



u 


Eo 


Ei 


E2 


E3 


E 4 


"7 


0.552275 


0.233202 


0.835423 


0.742113 


0.917366 


0.883988 


0.95098 


0.936032 


0.967737 


0.959912 


~T 


0.363492 


-0.033728 


0.796897 


0.670614 


0.903801 


0.861427 


0.944832 


0.92676 


0.964478 


0.955329 





-0.005928 


-0.384837 


0.743034 


0.572097 


0.886487 


0.832513 


0.937382 


0.91555 


0.960658 


0.967177 


7 


-0.83722 


-0.911674 


0.666155 


0.432269 


0.864211 


0.794284 


0.928328 


0.901662 


0.956183 


0.943688 








0.552275 


0.233202 


0.835423 


0.742113 


0.917366 


0.883988 


0.95098 


0.936032 


c 


+ 




+ 




+ 




+ 




+ 





For comparison, we give graphs of five low energy levels (Z = 180, j = 1/2) as functions of v, for ( = + (FIG. [5]i), 
for £ = — (FIG. [SJd). In addition, to give an idea of energy scale, we give a Table 4 of numerical values of some energy 
levels. 
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FIG. 4: ^-dependence of energy levels E„ (138, 1/2, f = ±) and Z-dependence of z/_ 



b=i/2 ' 



1.0- 

E/m 




C = + 
(a) 









. n=2 

n=1 






c=- 




(b) 





n n 

2 " 2 



FIG. 5: ^-dependence of energy levels (180, 1/2, C = ±) 



CONCLUDING REMARKS 



Here, we list the properties of all the radial s.a. Hamiltonians which are common to both nonsingular and singular 
regions. 

i) The spectrum of each s.a. Hamiltonian h (Z, j, £) is simple (nondegenerate) and contains a continuous part that 
is the set (— oo, —to] U [to, oo) and a discrete part {E n (Z,j, ()} located in the interval [— m, m); for a precise meaning 
of a nonnegative integer n, see the corresponding subsection. 

ii) The discrete spectrum is always accumulated at the point E — m, and the asymptotic form of the difference 
^nonrci _ ^ (Z,j, () — m as in — > oo is given by the well-known nonrelativistic formula 

^ onrcl = -mq 2 {2n 2 y 1 , 

the nonrelativistic spectrum does not depend on the extension parameter, j and ( (is degenerate in j and (). 

iii) Eigenfunctions of the discrete spectrum and generalized eigenfunctions of the continuous spectrum form a 
complete orthonormalized system in L 2 (R + ). 

As soon as all s.a. radial Hamiltonians h (Z,j, Q are fixed unambiguously, a corresponding total s.a. Dirac Hamil- 
tonian H (Z) is defined in a unique way. 

Because s.a. radial Hamiltonians h (Z, j, £) arc unique for Z < 118, the total Dirac Hamiltonian H (Z) with Z < 118 
is defined uniquely. For Z > 119, there is a family {H„ lt ... tU& (Z)} of possible total s.a. Dirac Hamiltonians. The 
family is parametrized by the parameters v± € [— n/2, n/2], —ir/2 ~ n/2, i = 1, A. The number A of the parameters 
is given by A = 2k(Z), where the integer k(Z) is given by k(Z) = (1/4 + Z 2 o?) x l 2 — 5, < S < 1. Any specific 
s.a. Dirac Hamiltonian H Vl ^„ tl/& (Z) corresponds to a certain prescription for a behavior of an electron at the origin. 
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The general theory thus describes all the possibilities that can be offered to a physicist for his choice. This choice 
is a completely physical problem. We believe that each s.a. Dirac Hamiltonian with superstrong Coulomb field can 
be understood through an appropriate regularization of the potential and a subsequent limit process of removing the 
regularization. We recall that a physical interest in the electronic structure of superheavy atoms was mainly motivated 
by a possible pair creation in the superstrong Coulomb field. Consideration of this effect in the framework of the 
most simplest model of a point-like nucleus was accepted to be impossible due to the conclusion (which is wrong as it 
is clear now) that this model is mathematically inconsistent [7] . We believe that a rehabilitation of the model allows 
returning to a consideration of the particle creation in this model providing considerable scope for analytical studies. 
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